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Abstract 

This  article  analyzes  the  problem  about  the  missile  overload  requirement  in  a  homing  terminal  guidance  under  various  en¬ 
gagement  scenarios.  An  augmented  proportional  navigation  guidance  (APNG)  model  is  introduced  on  the  basis  of  linear  kine¬ 
matics.  To  analyze  the  peak-to-peak  performance  of  the  terminal  guidance  system,  a  new  finite  time  L\  performance  measure  for 
a  linear  time- varying  (LTV)  continuous  system  is  proposed.  Then,  according  to  the  idea  of  the  adjoint  system,  a  novel  method  for 
computing  the  Lx  norm  of  a  linear  continuous  system  is  first  derived.  Within  the  finite  time  L\  framework,  the  quantitative  rela¬ 
tion  between  the  guidance  loop  dynamics  and  the  maximum  missile-target  maneuver  ratio  is  offered.  This  relation  is  expressed  in 
the  form  of  graphs  and  formulas  that  can  be  used  to  synthesize  some  of  the  major  subsystem  specifications  for  the  missile  guid¬ 
ance  system.  The  illustrative  examples  show  that  a  significant  performance  improvement  is  achieved  with  the  proposed  guidance 
loop  dynamics. 
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1.  Introduction 

In  the  missile  guidance  process,  the  missile  overload 
in  the  terminal  phase  is  critical  to  the  resulting  en¬ 
gagement  outcome.  In  practice,  missile  subsystems, 
such  as  autopilot  and  guidance  systems,  are  usually 
designed  separately.  After  combining  them  together, 
analysis  of  missile  overload  requirement  would  be¬ 
come  very  difficult  and  is  thus  a  severe  challenge  to 
control  engineers  because  of  their  interactive  effects 
and  uncertainties11"41. 

In  the  missile  terminal  guidance  process,  the  exter¬ 
nal  input  of  the  guidance  system  is  the  target  maneuver 
with  limited  acceleration  amplitude.  Moreover,  the 
major  source  of  the  missile  overload  is  the  target  ma¬ 
neuver.  In  engineering  practice,  to  design  nonsaturat¬ 
ing  guidance  systems,  the  guidance  system  designers 
need  to  know  the  missile  maximum  overload  require¬ 
ment  against  the  target  maneuver151.  Therefore,  the 
maximum  missile-target  maneuver  ratio  under  various 
engagement  scenarios  should  be  analyzed  since  the 
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maximum  acceleration  of  the  target  is  relatively  easier 
to  be  estimated  than  other  maneuver  information  of  the 
target.  In  Ref.[5],  the  maximum  missile-target  maneu¬ 
ver  ratio  is  obtained  under  some  assumptions  on  the 
guidance  loop  dynamics.  This  study  proposes  a  theo¬ 
retical  and  quantitative  approach  to  analyze  the  missile 
maximum  overload  requirement  for  the  worst-case 
target  maneuver  without  any  assumption. 

It  is  also  noted  that  there  are  some  salient  features  in 
the  terminal  guidance  analysis  quite  different  from  the 
conventional  performance  analysis.  First,  the  plant  is 
linear  time-varying  (LTV)  and  possibly  nonlinear 
while  in  most  performance  analyses16"81,  since  the 
nominal  plant  is  assumed  linear  time-invariant  (LTI), 
for  the  nonlinear  and  time- varying  effects  are  absorbed 
by  the  uncertainty  block.  Second,  missile  guidance  is 
such  a  complicated  control  problem  that  only  the  sys¬ 
tem  behaviors  in  a  finite  time  interval  should  be  con¬ 
sidered191.  It  means  that  these  special  physical  prob¬ 
lems  need  to  be  solved  within  a  new  framework. 

First  introduced  by  M.  Vidyasagar1101,  the  L\  ap¬ 
proach  formulated  the  problem  of  optimal  disturbance 
rejection  in  the  case  of  the  disturbance  which  gener¬ 
ated  as  an  output  of  a  stable  system  in  response  to  an 
input  which  was  assumed  to  be  either  of  unit  ampli¬ 
tude  or  arbitrary.  This  approach  is  of  considerable 
practical  importance  because  it  aims  to  minimize  the 
maximum  value  of  the  system  error.  In  the  classical 


control  system  performance  analysis,  such  as  L.,  norm 
minimization,  an  underlying  assumption  is  that  the 
system  external  inputs  have  their  power  in  finite  fre¬ 
quency  ranges  and  are  bounded  energy  signals.  From 
this  point  of  view,  the  L\  approach  serves  as  a  com¬ 
plement  to  the  theory  of  L„  norm  minimization. 

This  article  attempts  to  analyze  missile  overload  re¬ 
quirement  against  the  worst-case  target  maneuver 
based  on  the  guidance  loop  dynamics  over  a  finite  time 
interval  with  the  novel  finite  time  L\  approach.  First, 
the  L\  performance  measure  of  a  LTV  system  is  gener¬ 
alized  by  introducing  the  notion  of  finite  time  LTV 
system  behaviors.  Then,  an  effective  algorithm  of  the 
Zoo-induced  norm  for  a  LTV  continuous  system  is  first 
proposed  by  constructing  the  adjoint  system  that  asso¬ 
ciates  with  the  original  system  at  a  given  moment. 
Further,  this  algorithm  can  also  be  employed  to  the 
continuous  system  with  any  of  the  following  charac¬ 
teristics:  infinite  time,  finite  time,  LTI,  single-input- 
single-output  (SISO),  and  multi-input-multi-output 
(MIMO).  Based  on  the  LTV  terminal  guidance  model, 
the  quantitative  relation  between  the  maximum  mis¬ 
sile-target  maneuver  ratio  and  guidance  loop  dynamics, 
such  as  effective  navigation  ratio,  guidance  time  con¬ 
stant,  and  the  modified  coefficient  in  augmented  pro¬ 
portional  navigation  guidance  (APNG),  is  studied  with 
the  novel  finite  time  L\  approach.  In  addition,  the 
analysis  and  simulation  results  provide  the  feasible 
theoretical  basis  for  the  design  of  a  missile  terminal 
guidance  system. 

2.  Mathematical  Modeling 

Fig.l  shows  a  missile  terminal  guidance  loop  which 
typically  includes  the  dynamics  of  a  seeker,  a  guidance 
computer,  and  missile  dynamics1' IJ.  In  this  work,  the 
seeker  is  simply  viewed  as  a  pure  differentiator  which 
generates  the  line-of-sight  (LOS)  rate  of  the  target 
without  delay.  The  guidance  computer  is  viewed  as  an 
algebraic  function  that  computes  the  acceleration 
command  with  the  information  from  the  seeker.  Mis¬ 
sile/autopilot  dynamics  that  includes  vehicle  dynamics, 
autopilot,  and  actuator  is  viewed  as  a  simple  first-order 
or  second-order  system. 


Fig.l  Missile  terminal  guidance  system. 


The  general  formulation  of  a  3D  interception  prob¬ 
lem  is  rather  complicated.  However,  on  the  assumption 
that  the  lateral  and  longitudinal  maneuver  planes  are 
decoupled  by  means  of  roll  control,  the  3D  problem 
may  quite  reasonably  be  treated  as  an  equivalent  2D 


one.  Fig. 2  shows  the  kinematic  model  for  a  planar  in¬ 
terception  when  a  homing  missile  intercepts  a  maneu¬ 
vering  target,  where,  R  represents  the  relative  range 
between  the  target  and  missile  along  the  LOS,  q  is  the 
LOS  angle,  and  am  and  a{  are  the  missile  and  target 
accelerations  normal  to  the  LOS,  respectively.  The 
kinematic  relation  between  the  target  and  missile  mo¬ 
tion  can  be  described  by  the  nonlinear  equations: 

Rq  +  2Rq  =at-am  (1) 

R-Rq2  =alr-amT  (2) 

where  amr  and  atI  represent  the  missile  and  target  ac¬ 
celerations  along  the  LOS,  respectively.  If  the  missile 
approaches  the  target  with  a  constant  closing  speed  Vc, 
Eq.(2)  can  be  ignored  and  the  range  R  is  simply  ex¬ 
pressed  by  R=VC  fg0,  where  t„0  denotes  the  time  taken 
by  the  missile  to  intercept  the  target. 


Fig.2  Planar  missile  and  target  engagement  geometry. 

It  is  well  known  that  the  most  advanced  guidance 
law,  regarded  as  the  augmented  proportional  naviga¬ 
tion  guidance  law,  is  compounded  of  the  proportional 
navigation  guidance  law  and  the  target  acceleration 
compensation.  Then,  the  guidance  command  amc  is 
derived  from 

amc  =  NVcq  +  ^-at  (3) 

where  N  is  the  effective  navigation  ratio,  often  chosen 
to  be  3-5  in  practice,  and  K  is  the  modified  coefficient 
to  compensate  the  target  maneuver  normal  to  the  LOS. 

Missile/autopilot  dynamics  may  be  represented  by 

xm  =  Fmxm  +  Gmamc  (4) 

am=Hmxm  (5) 

or 

am=G(s)amc  (6) 

where  {Fm,  Gm,  Hm}  is  the  state-space  realization  of 
the  missile/autopilot  dynamics,  and  G(s)  denotes  the 
transfer  function. 

For  simplicity,  the  missile/autopilot  dynamics  is 
viewed  as  a  first-order  system,  but  it  should  be  pointed 
out  that  any  linear  high  order  missile/autopilot  system 
is  also  allowed  within  the  finite  time  L\  framework.  By 
combining  Eq.(l)  with  Eqs.(3)-(5),  the  APNG  loop 
dynamics  can  be  expressed  by  the  following  linear 


No.4 

Ji 

Denggao  et  ah 

time-varying  differential  equations: 

2  -1  1 

i 

r  4 1 

K-t  ^(*,-0  r*i 

NVC  -1  LflmJ 

K 

~Y~  Y  \ 

IT 

z  =  [0 

€1 

where  T  is  the  time  constant  of  the  missile/autopilot, 
which  is  equivalent  to  a  state-space  realization  with 
Fm=  -l IT,  Gm=\/T,  and  Hm=  1,  and  te  is  the  total  flight 
time  of  the  engagement.  Note  that  te  is  now  a  constant. 
The  value  of  (te-t)  is  the  time  to  go  until  the  end  of  the 
flight. 

From  Eq.(7),  it  is  noted  that  the  missile/target  en¬ 
gagement  in  the  terminal  phase  can  be  regarded  as  an 
LTV  system  in  the  finite  time  interval  [0,  te).  In  this 
article,  the  output  of  interest  is  the  missile  lateral  ac¬ 
celeration  throughout  the  flight. 


3.  Finite  time  L{  Performance  Analysis  of  LTV 
System 


This  section  will  present  the  finite  time  L\  perform¬ 
ance  measure  and  describe  the  way  it  is  calculated  for 
a  given  LTV  continuous  system. 


3.1.  Performance  measure 


In  this  subsection,  the  performance  measure,  the  Loo- 
induced  norm,  is  defined  for  an  LTV  continuous  sys¬ 
tem  with  a  special  emphasis  on  finite  horizon  behav¬ 
iors.  Various  definitions  for  LTV  systems  presented  in 
this  article  are  borrowed  from  Ref.  [12], 

Eq.(7)  can  easily  be  converted  into  the  canonical 
LTV  form  as  follows: 

£ .  { m  =  At)x(t)  +  B{t)w(t),  *(0)  =  0 

'!*(*)= C(#M0 

where  jc(f)eR"  is  the  system  state  vector,  w( f)  e  Z.2[0,?r] 
the  exogenous  disturbance  signal,  jc(0)  the  initial  state, 
j(i)eR"cR"  the  state  combination  (objective  function 
signal)  to  be  attenuated,  and  A(t),  B(t),  C(t)  the 
bounded  continuous  matrices. 

To  formalize  the  definition  of  the  finite  time  L\  per¬ 
formance  measure,  Definition  1  establishes  the  nota¬ 
tion  for  La,  signal  norms  over  the  finite  horizon  used 
throughout  the  article. 

Definition  1  For  a  finite  constant  t{>  0,  the  fi¬ 
nite-horizon  infinity  norm  of  a  signal  / :  R  —  R+  is 
defined  as 

(9) 

where  ||«||  denotes  the  Euclidean  L2  vector  norm.  The 
linear  space  L40,  L]  is  expressed  by 
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/-,[0,/r]:={/ :R*  >R|/,J0;1<x}  (10) 

The  subset  {/ :  R+  ->  R+  <  l}  c=  Lx  [0,tf  ] 

is  denoted  by  SLK[0,tf], 

For  a  fixed  time  t{  >  0,  the  closed- loop  Lx  perform¬ 
ance  measure  of  the  system  Eq.(8)  is  defined  as 

(U) 

Here,  under  the  fixed  initial  condition,  jc(0)  =  0,  Eq.(ll) 
defines  a  mapping  from  bounded  amplitude  inputs 
w(t)&BLx[0,  tf]  to  bounded  amplitude  penalty  outputs 
z(t)eBL.y[0,  tf]  and  a  relevant  performance  criterion  is 
the  peak-to-peak  or  Loo-induced  norm  of  this  mapping. 
In  the  LTI  case,  Eq.(l  1)  can  be  derived  from 

J(Z,tf)  =  \l\h(t)\h  (12) 

where  h(t)  is  the  known  unit  impulse  response  of  the 
LTI  system1131.  However,  in  the  LTV  case,  Eq.(12)  be¬ 
comes  complicated  as  follows: 

sup  D#,:4lr  (13) 

<W0' 

where  h(t,f)  is  the  impulse  response  of  the  LTV  system. 
Since  h(t,r)  has  no  general  analytic  expression,  Eq.(13) 
can  not  be  used  to  compute  the  Loo-induced  norm  for 
LTV  system  in  practice.  To  the  best  of  authors’ 
knowledge,  there  is  no  mention  in  the  literature  of  the 
computation  method  of  the  Loo-induced  norm  for  a 
LTV  system  at  present. 

3.2.  Computation  method  of  L^-induced  norm  for  LTV 


Let  ta  (0  <  t0<  tf)  be  the  moment  when  the  output  of 
the  system  is  of  particular  interest.  The  adjoint  re¬ 
sponse  of  system  Eq.(8)  at  t0  is  defined  as 

±x*A%)  =  AT(t0-tg)x*%s) 

d  h 

z°A%)  =  B\t0-t%)x°A\t%) 

jcadj(0)  =  C'(t0) 
where  AJ,  BJ  and  CT  denote  the  transpose  of  matrices 
A,  B  and  C,  respectively;  Jtadj  is  the  state  adjoint  re¬ 
sponse;  zadj  the  output  adjoint  response,  t„  the  time 
variable  of  the  system  Eq.(14).  System  Eq.(14)  can 
also  be  regarded  as  the  auxiliary  system  used  to  compute 
the  Loo-induced  norm  of  the  linear  continuous  system. 

With  the  aid  of  system  Eq.(14),  the  following  theo¬ 
rem  can  be  obtained. 

Theorem  1  The  output  of  system  Eq.(8)  at  the  given 
moment  t0  is 

z(fo)  =  J0'”(2adj(^o-^))TM^)d^  (15) 

Proof  The  output  of  system  Eq.(8)  at  the  given 
moment  t0  can  be  written  into 
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z(t0)  =  j‘°  C(tomo,r)B(r  )n>(z)dz  (16) 


where  <£( t0,z )  is  the  transition  matrix  of  system  Eq.(8). 
It  is  noted  that  the  transition  matrix  by  its  definition 

satisfies  the  following  equations: 

f  4>(t,z)  =  A(t)<Ht,z) 
ot 

(17) 

-^-0(t,z)  =  -0(t,z)A(z) 
dz 

(18) 

4>(t,t)  =  I 

(19) 

0(t,z)0(z,Z)  =  0(t,Z) 

(20) 

where  /  is  the  identity  matrix  with  the  appropriate  di¬ 
mension.  According  to  Eqs.(17)-(19),  it  can  be  proved 
that  the  output  and  the  state  of  the  adjoint  system 
Eq.(14)  are  expressed  by 

z*di(z)  =  BT(t0-z)<PT(t0,t0-z)CT(t0)  (21) 

and 

x°di(z)  =  <PT(t0,t0-z)CT(t0 )  (22) 

respectively,  where  re[0,  f0],  By  combining  Eq.(16) 
and  Eq.(21),  Eq.(15)  can  be  obtained.  Thus,  the  proof 
is  completed. 

Based  on  Theorem  1,  this  article  proposes  the  main 
tool  to  compute  the  Too-induced  norm  of  the  LTV  con¬ 
tinuous  system  as  follows. 

Theorem  2  The  closed-loop  L\  performance  meas¬ 
ure  J(E,  if)  of  system  Eq.(8)  must  satisfy  the  following 
equation: 

j(Z,t{)=  sup  f  ,o||zadj(t0^f||dz-  (23) 
<„do,<f]JO  11  W-' 

where  1*^  is  the  column-sum  norm  of  a  matrix. 

The  proof  of  Theorem  2  is  simple  in  the  SISO  case. 
Combining  Eq.(ll)  with  Eq.(15)  results  in  the  Eq.(23). 
It  is  also  noted  that  when  w(z)  in  Eq.(15)  satisfies 
Eq.(24),  the  finite  time  infinite  norm  of  the  system 
output  z(t)  reaches  the  maximum. 

w(t)  =  sgn(zadj(/0  —  z))  (24) 

where  sgn(  • )  is  a  signum  function  that  extracts  the 
sign  of  a  real  number  and  re[0,  /<,]• 

It  is  important  to  point  out  that,  a  significant  advan¬ 
tage  over  Eq.(13)  lies  in  that  the  value  of  jadj( t0  -  z)  in 
Eq.(23)  can  be  easily  acquired  through  the  auxiliary 
system  Eq.(14).  Moreover,  Eq.(23)  is  suitable  for  the 
linear  continuous  system  in  such  cases  as  finite  hori¬ 
zon,  infinite  horizon,  LTI,  LTV,  SISO,  and  MIMO. 
Also,  in  the  finite  horizon  case,  the  system  matrix  A(t) 
in  Eq.(8)  may  not  be  stable.  Therefore,  the  proposed 
finite  time  L\  approach  is  quite  fit  to  analyze  the  mis¬ 
sile  terminal  guidance  system  described  by  Eq.(7) 
since  it  turns  unstable  when  the  flight  time  t  ap¬ 
proaches  the  final  intercepting  time  tf. 

For  the  LTI  and  SISO  system,  Eq.(23)  can  be  sim¬ 
plified  into 
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y(27,ff)  =  j^j|#j(ff-r||t  (25) 

Eq.(25)  provides  another  method  for  computing  the 
L\  norm  of  a  LTI  system  as  well  as  for  solving  Eq.(12). 

4.  Missile  Overload  Requirement  Analysis 


The  concept  of  the  finite  time  Lx  approach  proposed 
in  Section  3  is  of  practical  use  to  analyze  the  missile 
overload  requirement  in  tackling  the  worst-case  target 
maneuver  since  the  guidance  loop  operates  in  a  finite 
time  interval.  Consequently,  it  is  necessary  to  analyze 
the  relation  between  the  guidance  loop  dynamics  and 
the  upper  bound  /./  on  the  required  missile-target  ma¬ 
neuver  ratio  defined  as 


B~ 


sup  |am(0| 


sup  |at(0| 

rqo*] 


(26) 


In  Ref.  [4],  let  ZG(s)  range  from  0°  to  180°  in 
Eq.(6),  then  the  maximum  missile-target  maneuver 
ratio  //  is  derived  from 

ft  =  V  ,VN  >  2  (27) 

H  N—2  y  J 

However,  to  alleviate  the  negative  effects  of  a  high 
frequency  noise,  the  phase  of  G(s)  in  high  frequencies 
should  usually  be  required  to  drop  rapidly  under  zero. 
In  this  respect,  contrary  to  the  existing  literature,  the 
proposed  finite  time  L\  approach  can  calculate  the 
maximum  missile-target  maneuver  ratio  without  any 
assumption  on  G(s)  and  the  detailed  form  of  the  worst 
target  maneuver  can  also  be  obtained  abiding  by 
Eq.(24). 

According  to  Eq.(l  1)  and  Theorem  2,  in  order  to  get 
the  maximum  missile-target  maneuver  ratio  without 
loss  of  generality,  only  the  normalized  unit  target  ac¬ 
celeration  needs  to  be  taken  into  account. 

The  typical  missile  terminal  guidance  loop  dynamics 
is  described  by  Eq.(7)  with  the  ensuing  coefficients: 
t{=  6  s,  Vc=\  500  m/s,  K=  0,  T=  0.3  s,  N=  4. 

In  Fig.3,  assuming  that  the  guidance  dead  zone  is 
300  m,  the  worst-case  target  maneuver  form  is  shown 
over  the  finite  time  interval  [0,  5.8].  The  vertical  axis 
denotes  the  normalized  acceleration  of  the  target.  The 
missile  maximum  overload  occurs  at  5.8  s  when  the 
target  performs  the  worst  maneuver  as  shown  in  Fig.3. 

Fig.4  shows  the  quantitative  relation  between  the 
maximum  missile-target  maneuver  ratio  and  the  effec¬ 
tive  navigation  ratio  when  T=  0.3  s.  To  avoid  the  lar¬ 
ger  missile  overload  requirement,  N  is  expected  to  be 
in  the  interval  [3.0,  4.5].  Compared  to  the  results  based 
on  the  conventional  linear-quadratic  performance  in¬ 
dex  with  terminal  constraints,  Fig.4  shows  not  only  the 
optimal  interval  for  the  effective  navigation  ratio,  but 
also  the  quantitative  maximum  missile-target  maneu¬ 
ver  ratio  as  well.  Since  the  maximum  target  maneuver 
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ability  is  often  known  in  advance,  the  missile  maxi¬ 
mum  overload  requirement  can  be  obtained  with  the 
finite  time  Lx  approach. 


Flight  time/s 


Fig.3  History  of  the  normalized  acceleration  of  the  worst- 
case  target  maneuver. 


Effective  navigation  ratio 


Fig.4  Maximum  missile-target  maneuver  ratio  with  differ¬ 
ent  effective  navigation  ratios. 

Fig. 5  shows  that  the  guidance  time  constant  T  is  the 
main  contributor  to  the  missile  overload,  and  the 
maximum  missile-target  maneuver  ratio  is  proportional 
to  the  guidance  time  constant.  When  T  approaches 
zero,  the  maximum  missile-target  maneuver  ratio 
equals  2,  which  accords  with  the  result  of  Eq.(27) 
when  N=  4.  This  indicates  that  reducing  the  response 
time  of  the  missile/autopilot  system  not  only  delays  the 
destabilizing  moment  of  the  terminal  guidance  system 
but  also  decreases  the  missile  overload  requirement. 


Guidance  lime  constant/s 

Fig.5  Maximum  missile-target  maneuver  ratio  with  differ¬ 
ent  guidance  time  constants. 


Fig. 6  clearly  depicts  that  APNG  is  superior  to  the 
conventional  proportional  navigation  guidance  (K=  0) 
in  that  an  appropriate  modified  coefficient  to  compen¬ 
sate  the  target  maneuver  can  be  chosen  if  the  real-time 
estimated  target  acceleration  is  available,  which  leads 
to  an  enormous  reduction  of  missile  overload  require¬ 
ments.  Fig. 6  shows  the  optimal  modified  coefficient  in 
APNG  is  1.5. 


Modified  coefficient  of  APNG 


Fig.6  Maximum  missile-target  maneuver  ratio  with  differ¬ 
ent  modified  coefficient  of  APNG. 

Fig. 7  compares  three  cases  of  the  guidance  loop 
dynamics.  In  Case  1,  K,  N  and  T  equal  0,  4,  0.3  s,  re¬ 
spectively.  In  Case  2,  all  parameters  remain  unchanged 
excepting  the  guidance  time  constant  T  which  is  0.02  s. 
In  Case  3,  the  better  modified  coefficients  are  chosen 
according  to  Figs.4-6,  where  K  is  1.5.  It  should  be 
stressed  that  the  history  of  the  target  maneuver  varies 
with  flight  time.  When  a  proper  guidance  loop  dynam¬ 
ics  is  designed,  the  maximum  missile-target  maneuver 
ratio  can  be  decreased  to  1.3. 


Flight  time/s 


Fig.7  History  of  the  maximum  missile-target  maneuver  ratio. 

From  Figs.3-7,  it  is  noted  that  the  proper  assessment 
of  /./  together  with  a  proper  design  of  the  guidance  loop 
dynamics  may  lead  to  a  satisfactory  nonsaturating 
guidance  system. 

5.  Conclusions 

A  finite  time  Lx  performance  measure  is  proposed 
for  a  LTV  system,  based  on  which  the  performance  of 


terminal  guidance  process  with  the  various  guidance 
loop  dynamics  is  analyzed.  The  missile  overload 
requirement  for  the  terminal  guidance  is  obtained  and 
some  design  guidelines  are  given.  It  is  emphasized  that, 
within  the  proposed  novel  L\  framework,  the  LTV  con¬ 
tinuous  system  is  not  required  to  be  stable  which  is 
suitable  for  analysis  of  terminal  guidance  system. 
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